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ABSTRACT 

Context. High-energy gamma-ray emission from the Galactic plane above ~ 100 MeV is composed of three main contributions: diffuse 
emission from cosmic ray interactions in the interstellar medium, emission from extended sources, such as supernova remnants and 
pulsar wind nebulae, and emission from isolated compact source populations. 

Aims. The diffuse emission and emission from the extended sources provide the dominant contribution to the flux almost everywhere 
in the inner Galaxy, preventing the detection of isolated compact sources. In spite of this difficulty, compact sources in the Galactic 
plane can be singled out based on the variability properties of their y-ray emission. Our aim is to find sources in the Fermi data that 
show long-term variability. 

Methods. We performed a systematic study of the emission variability from the Galactic plane, by constructing the variability maps. 
Results. We find that emission from several directions along the Galactic plane is significantly variable on a time scale of months. 
These directions include, in addition to known variable Galactic sources and background blazars, the Galactic ridge region at positive 
Galactic longitudes and several regions containing young pulsars. We argue that variability on the time scale of months may be 
common to pulsars, originating from the inner parts of pulsar wind nebulae, similarly to what is observed in the Crab pulsar. 

Key words. Survey s,Gamma rays: general, Methods: data analysis 



1. Introduction 

Most of the y-ray emission from the Galactic plane in the en- 
ergy band above 0.1 GeV is produ ced by cosmic ray interac- 
tions with interstellar matter (Ficht eTet al.||1975] |Hunter et al.| 
|1 9971 1 Abdo et al.|2009a| ). These interactions, on a time scale of 
~ 10 v yr and on distance scales of hundreds of parsecs, result 
in the production of bright large-scale diffuse emission from the 
entire Galaxy. Isolated Galactic y-ray sources, such as pulsars, 
pulsar wind nebulae, supernova remnants, and y-ray-loud binary 
systems are superimposed on this large-scale diffuse emission. 
The strong diffuse background emission provides the main ob- 
stacle for detecting individual isolated sources in the Galactic 
plane. Inhomogeneities of the diffuse emission induced by the 
inhomogeneity of matter distribution in the interstellar medium 
lead to local brightness enhancements of the y-ray emission from 
the Galactic plane, which might be misinterpreted as isolated 
sources. Some of the sources at low Galactic latitude, reported 



in the tw o-year Fermi catalog, might have such a nature ( [Nolan 
|et aL|2Qll| ). 

y-ray-loud binaries provide a very moderate contribution to 
the overall y-ray emission from the Galaxy. Up to now, only a 
few binary systems were found to produce y-rays with energies 
above 0.1 GeV. Six out of nine known y-ray-loud binaries, PSR 
B 125 9-63 ( |Aharonian eTaLl|2005]), LSI +61 303 ( [Albert et al.l 
2006]), LS 5039 ( [Aharonian et al.|2006] ), Cyg X-3 ( [Abdo et"aL] 
2009b] ), Cyg X-l ( [Albert et al.||2007| ), and HESS J0632+057 



( Casares et al.|2012"| ) are systems composed of a compact object 
(a neutron star or a black hole) orbiting a m assive companion 
star. One system, Eta Car (Abd o et al.|20 10b), is a massive star 
binary with Wolf-Rayet star, one system, V407 Cyg ([Abdo et aTf 



2010c]), is a symbiotic binary containing a white dwarf accreting 
matter from the main sequence star. The nature of the recently 
disc overed source 1FGL J1 01 8.6-5856 is not well constrained 
yet ( [Ackermann et al.|2012 ). 

Indeed, only a very small fraction of the known high-mass 
binaries with a compact object orbiting a massive star (most of 
them are discovered as high-mass X-ray binaries, HMXRB) turn 
out to be "y-ray-loud". It is not clear at the moment if these 
are only few detected y-ray-loud binary systems because only 
a small fraction of the HMXRB s are able to produce y-rays 
or because most of the HMXRB s are emitting y-rays, but their 
typical flux levels are just below the sensitivity of Fermi. It is 
also not clear which properties of the HMXRB are responsible 
for the enhanced y-ray emission from the four known y-ray- 
loud HMXRB s. The distinguishing feature of the y-ray-loud bi- 
nary systems is variability. Other known types of Galactic y-ray 
sources, including the diffuse emission and supernova remnants, 
are not variable on day-to-month time scales. This feature could 
reveal the presence of y-ray -loud binaries even at the top of 
much stronger diffuse emission. 

Among pulsars, only the Crab pulsar, which is the youngest 
radio pulsar in the Galaxy, is known to produce y-ray flares on 
time scale of hours-to-days ( Abdo et al. 201 1| ). The variable 
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emission is produced in the inner part of the pulsar wind nebula 
(PWN) associated to the Crab pulsar. It is not clear if the flaring 
activity of the nebula is a peculiar feature of the Crab pulsar or if 
it is a generic property of the young-pulsar-powered systems. A 
systematic study of variability properties of other known young 
pulsar systems in the Galaxy is needed to clarify this question. 

Fermi/LAT is a pair-conversion gamma-ray detector operat- 
ing between 20 MeV and 300 GeV. The LAT has a wide field of 
view of sr at 1 GeV, and observes the entire sky every two orbits 
( [Atwood et al.|2009] ). Therefore it is a perfect telescope for the 
systematical study of the variability of the GeV sky on a scale of 
hours-to-days. 



2. Search for variability of y-ray emission in the 
Galactic plane 

A simple way to verify if y-ray emission from a given direction 
on the sky is variable or constant is to analyze the lightcurve and 
check if it is consistent with a constant flux. Several measures of 
(in)consistency with the constant flux can be considered. 

Consider the flux F(l, b, U) = F;(Z, b) in the given sky di- 
rection (/,/?) (in Galactic coordinates) binned in N time bins 
1 < i < N. The flux is consistent with constant if the reduced 
X 2 of the fit of the lightcurve with the constant 



1 A (Fj-F(l,b)) 2 



(1) 



is close to one (here <x; are the measurement errors in each time 
bin and F(l, b) is the average flux in the (/, b) direction). If we 
consider the dependence of x 2 on (/, b) we expect to find x 2 - 1 
in all directions except for the directions in which there are 
variable isolated y-ray sources. All non- variable y-ray sources, 
such as supernova remnants and the large-scale diffuse Galactic 
emission should not appear in the x 2 (h b) map. Instead, the 
sources that are characterized by a variability amplitude larger 
than ~ o~i/F(l, b) should be visible as excess in the^ 2 map. 

The scatter of the source flux around its average value, char- 
acterized by the x 2 of the fit of the lightcurve with a constant, 
provides a good measure of variability especially for sources that 
are systematically variable throughout the observation period. 
On the other hand, sources that produce on average one flare 
with a duration much shorter than the span of the observation 
period are gradually "washed out" from the x 2 (h b) map as the 
overall exposure time increases. Indeed, the contribution of the 
time bins containing the flare to the overall^ 2 given by Eq. (]]) 
decreases with increasing number of bins. Accordingly, single- 
flare sources might disappear from the x 2 map with large expo- 
sure. A better measure of variability for the single-flare sources 
is the maximum deviation of the flux from the average value 
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Similarly to the^ 2 measure, excess of cr max beyond the expected 
statistical fluctuations indicates the presence of a flaring source 
in a given sky direction. In contrast to the x 2 measure, the cr max 
measure does not decrease with increasing exposure. 

A systematic view on the entire population of variable 
sources in the Galaxy (both systematically variable and flar- 
ing) can be obtained by identifying all significant excesses in 
the x 2 (L b) and cr max (l, b) maps constructed based on the y-ray 



lightcurves in different energy bands and different energy scales. 
In the following sections we investigate these sets of "variability 
maps" of the y-ray emission from the Galaxy. 

Collecting the photons within a PSF (point spread function) 
can lead to false-positive or misaligned detections of variability. 
Consider the case of N strongly variable sources with maximum 
flux deviations from constant level at <x;, i = I..N. These sources 
can cause an excess of variability at some region, at which 
the contributions from the PSF of all these sources are non- 
negligible. Namely, if P(r) is the PSF (the fraction of source's 
photons beyond radius r) of Fermi/LAT and at is the angle un- 
der which the region is visible from i-th variable source position, 
it can be shown that the significance of false-variability of this 
region can be estimated as 



V false 



(3) 



where r t is the angular distance from i-th variable source to the 
selected region. 

Consider now the case of a flaring weak source in the vicinity 
of bright non- variable source. In this case the non-variable bright 
source will contribute "background" photons at the position of 
the weak source. This background reduces the significance of 
the flare in photons collected within the PSF exactly around the 
position of the flare. However, the number of "background" pho- 
tons will decrease with increasing distance from the bright non- 
variable source. This can shift the observed maximum of the 
variability from the flare position along the line connecting a 
bright and a flaring source in the direction farther from the non- 
variable source. 

Possible false detections of variability caused by the influ- 
ence of nearby sources are best addressed using variability maps: 
the imaging of the variability properties of emission helps to 
spot these possible false detections of variability of non-variable 
sources in the presence of known nearby variable sources. 

3. Observations and data analysis 

In our analysis we used all publicly available LAT data for the 
period from August, 2008 to October, 2011. We processed the 
data using the Fermi Science Tools software package, version 
v9r23pl [] We filtered the event lists using gtselect tool and left 
only "source" events (evclass=3) at zenith angles 6 Z < 100. 

To produce the variability maps, we used the "aperture pho- 
tometry" method. We binned the events from circles of radius 2° 
around each reference position on the sky, using the gtbin tool. 
The exposure in each time bin and each sky direction was cal- 
culated using the gtexposure tool and the instrument response 
function set for the eventclass P7S0URCE_V6. 



4. Results 

4.1. Variability of emission from the Inner part of the Galactic 
disk 

The lightcurves with month-long time bins for the 3.1 yr expo- 
sure of LAT have 37 time bins. One can estimate the probability 
for cr max in any individual lightcurve to exceed a certain thresh- 
old value cr thr . For example, the chance probability for cr max to 
exceed the threshold of <x^ r = 3.2 is 5%, while only 0.2% of the 
non-variable source lightcurves would occasionally have cr max 



http : //f ermi . gsf c . nasa . gov/ssc/data/analysis/ 
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Fig. 1. LAT count map and variability contours {cr max values) 
in the energy band above 300 MeV in the directions of the sky 
regions containing the LSI +61 303 y-ray-loud binary (left) and 
the Cyg X-3 and V407 Cyg y-ray-loud binaries (right). Contour 
levels start from 3<x with steps of 0.5<x 

in excess of cr thr = 4. Producing the variability map of the en- 
tire Galactic plane, we have scanned over all different directions 
along the 360° of the plane projection on the sky. The PSF of 
LAT at energies above 300 MeV has a half- width (68% contain- 
ment radius) of 2°. This implies that ^ 360° /4° = 90 indepen- 
dent lightcurves were probed. If the flux from the entire Galactic 
plane is constant, one expects to find less than one lightcurve (on 
average) in which cr max is in excess of cr t h r = 4. Consequently, 
all the directions in which variability in excess of cr thr = 4 is 
found should be considered as directions from which real vari- 
able emission is detected. 

Fig. [T] shows the LAT count map in the energy range above 
300 MeV with the superimposed contours showing the cr max val- 
ues found in the lightcurves binned into the month-long time 
bins for the sky regions containing the known variable Galactic 
y-ray sources LSI +61 303, Cyg X-3 and V407 Cyg. All three 
sources are y-ray-loud binaries with quite different types of vari- 
ability. LSI +61 303 is a persistent source, variable on 26.4 d 
orbital period and on a longer superorbital period of 4.6 yr. Cyg 
X-3 is a recurrent flaring source that already exhibited several 
bright flares over the 3.5 yr period of Fermi/LAT exposure. V407 
Cyg is a symbiotic binary that produced only one flare during the 
LAT exposure. One can see that in spite of the large difference 
in the variability properties of the sources, all three sources are 
clearly detected as strong excesses in the variability map, shown 
by the green contours. The other excesses above and below the 
Galactic plane are produced by known blazars. 

Fig. [2] shows the LAT count map with superimposed vari- 
ability contours for the region around the Galactic Center. The 
excess variability is not centered at the Galactic Center itself (the 
bright source in the center of the image). Instead it traces the part 
of the Galactic plane at positive Galactic longitude 0° < / < 7° 
dominated by diffuse emission, with no clearly isolated bright 
sources. The end part of the variable region, with the brightest 
excess at the level of cr max 4.5, is centered at the bright isolated 
source identified with supernova remnant W28. Appearance of 
excess variability from the region of diffuse emission and from 
a supernova remnant is surprising. 

Another region with excess variability is found in the direc- 
tion of the Kookaburra pulsar wind nebula complex, shown in 
Fig. [3] This sky region contains three young pulsars, PSR J1418- 
6068, PSR J1410-6132 and PSR J1413-6205, which emit pulsed 
y-ray emission in the 0.1-10 GeV band. Similarly to the regions 
of the Galactic ridge and of LSI +61 303, Cyg X-3 / V407 Cyg, 
the variable sources above and below the Galactic plane, visible 
on the map, are background blazars. The source to the right of 
the Kookaburra complex, PSR B 1259-63, is not visible in the 
three-year count map (shown by the color scale), but appears 




Fig. 2. LAT count map and variability contours in the direction 
of Galactic ridge and W28 supernova remnant. Notations are the 
same as in Fig.[T] 




15 29 43 57 71 84 98 112 



Fig. 3. LAT count map and variability contours in the direction 
of sky region containing pulsars PSR J141 8-6068, PSR J1410- 
6132 and PSR J1413-6205. Notations are the same as in Fig. [2] 

as a very bright excess in the variability map. This is not sur- 
prising, because this source is a y-ray-loud binary with a radio 
pulsar on a highly eccentric orbit with a 3.5 yr orbital period. It 
exhibits periods of activity only close to the periastron passage. 
The last periastron passage happened in December 2010. The 
source produced a bright one-month-long flare in January 2011 
and disappeared since then. 

Excess of variability on a scale of months at the level cr max 
15 is also detected from the direction of the Crab pulsar. 

5. Discussion 

Emission variability from the Galactic plane is caused by com- 
pact y-ray sources. The only firmly established class of generi- 
cally variable Galactic sources is that of y-ray-loud binaries. In 
this respect, it is not surprising that the highest excess variability 
is found in the directions of known y-ray-loud binary systems: 
LSI +61 303, V407 Cyg, Cyg X-3, and PSR B 1259-63. 

The recent discovery of y-ray flaring activity of the Crab pul- 
sar shows that pulsars (or the inner compact parts of the pulsar 
wind nebulae, PWNe) should also be considered as a possible 
variable Galactic source population. However, it is difficult to 
draw definitive conclusions about the variability of the entire 
pulsar population based on the study of a single source, espe- 
cially taking into account the fact that the Crab pulsar could be 
considered as a special case because it is the youngest known 
y-ray-loud pulsar. Our study shows, however, that long (month) 
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time scale variability may indeed be a generic feature of pul- 
sars/PWN systems, as discussed below. 

In particular, significant variability is found in the direction 
of the region of the Kookaburra nebula, which contains several 
young pulsars. Unfortunately, the point spread function of LAT 
in the energy band in which the variability is detected is too large 
to tell if one of the pulsars produces variable emission on a scale 
of months, or if the variability is the result of the variability of 
several pulsars. It is also not clear why pulsars in the Kookaburra 
region would be variable while other young y-ray-loud pulsars 
in the Galactic plane do not reveal any excess variability a time 
scale of months. Additional study is needed to clarify the origin 
of y-ray emission from young pulsars, including the youngest 
Crab pulsar. 

Another region that exhibits excess variability is the Galactic 
ridge at positive Galactic longitudes 0° < / < 7°. The strongest 
excess of variability is centered at the supernova remnant W28. 
Clearly, the variability in this sky direction cannot be produced 
by the extended emission from parts of W28 itself, which have 
extension of about 0.5° ( |Abdo et~aL||2010a| ), i.e. ~ 20 pc at 
~ 2 kpc. The variable region contains several young pulsars, 
including PSR J1747-2809 and PSR J1746-2850, close to the 
Galactic Center with the ages T = 5.3 kyr and T = 12.7 kyr, just 
slightly older than those of the Crab pulsar. Closer to the W28 
remnant there is PSR B 1757-24, a pulsar of age T = 15.5 kyr. 
Finally, W28 its elf was su pposed to be associated with the pulsar 
PSR 1758-23 ( [Frail et al.|1993) . It is possible that the variability 
is produced by one or several pulsars in this region of the Galaxy. 
Otherwise, the source of variability could be y-ray-loud binaries 
that are below the sensitivity of LAT. 

We also remark that the position of this variable region is 
consistent with the positions of sources that are variable in hard 
X-ray band (see e.g. [Tblezhinsk y'etaLl ([20 10)). 

The threshold value cr thr for significant detection of variabil- 
ity is different for the search of variable emission from an arbi- 
trary direction along the Galactic plane and for the search of vari- 
ability of emission from known isolated y-ray source classes. If 
we are interested in the variability properties of known sources, 
the threshold cr thr can be decreased, because in this case there is 
no trial factor associated to the scan of the entire Galactic plane. 
The probability to find an excess variability at the level above 
3.5cr is just about 1% for each individual source. 

Excess variability at a level higher than 3.5<x is found in two 
more directions on the sky, which coincide with the positions 
of the known sources, PSR J1826-1256 and PSR Jl 119-6127. 
The LAT count maps with the superimposed variability contours 
in these directions are shown in Fig. [4] PSR J1826-1256 is sit- 
uated close to a known y-ray-loud binary LS 5039, which ex- 
hibits emission variable on the 4.6 d orbital period. Surprisingly, 
the emission from LS 5039 appears to be very stable on a time 
scale of months, so that no excess variability associated with 
LS 5039 is found. No variability is detected either from the two 
other binary systems, Eta Car and 2FGL J1019. 0-5856 (see right 
panel of Fig. |4j. Evidence for excess variability at the locations 
of these pulsars, if combined with the detection of variability 
from the direction of Kookaburra region that contains several 
young pulsars, shows that long time scale variability may be a 
generic property of young pulsars, and not a peculiar feature of 
the Crab pulsar alone. 

To summarize, we have studied variability of y-ray emission 
on a time scale of months from the Galactic plane in the energy 
band above 300 MeV in the Fermi/LAT data and found several 
directions with significantly variable emission. We found that 
strong variability excess is detected in the directions of known y- 




Fig. 4. Variability contours superimposed on the LAT count map 
for the sky region of HESS Jl 825- 137. Notations are the same 
as in Fig. [I] 

ray-loud binary systems. Remarkably, there are no signatures of 
weaker variability, which could be clearly associated to other bi- 
nary systems with black holes and/or neutron stars. At the same 
time, our study revealed weaker variability excess in the direc- 
tions of pulsars. This indicates that long variability on a scale 
of months, as recently discovered in the Crab pulsar, may be 
a generic feature of the y-ray emission from the young pulsar 
population, so that pulsar/PWN systems should be considered as 
variable Galactic y-ray sources. 
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